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1 . INTRODUCTION 

Certainly you have met the number tt frequently in 
*yotir study of mathematics. You have learned various - 
approximations^ i-t. • Perhaps the most familiar are 
^ 22/7 and 3 .,14. Through the centuries there have' been 
many other approximations of tt . In 1967 a computer 
program calculated Tr^to 500,000 decimal places/ But no 
one has ever f found an "exact" value of tt , and no one 
ever will, for a very good reason. £he title of this 
unit states that reason, and in this unit we shall 
* prove the statement-. 

l^lV Rational and Irrational Numbers 

r 

Recall that, a real number is said to be rational 
if#ami pnly' if i£ can be expressed in the form p/q " 
(i.$. , as- a ratio) where p and q are integers. A real 
number which is not rational is said to be irrational. 
On the face of it, it is not clear that there are any. 
^irrational nyfcbers. In fact, the ancient Greeks believed 
tha.t all numbers were rational. We can- hardly blame therif, 
since the Greeks thought geometrically, and geometric 
"common sense" seems to "confirm their belief. After all, 
we can draw a straight -line "(in modern terms, think of the 
^ x-axis) and mark off upon it equally spaced points corre- 
sponding to the integers. 'Then, by a straightforward tech- 
nique of Euclidean geometry, we can divide each unit inter- 
val on the line into q equal intervals, where q is any 
positive integer we wish. If we, do this, for example, to 
the interval froto 0 ttT^we get points whose distances from 

0 are \> f> : • •> If we do this for all possible q 

and for all unit interv als, it is hard for us (and was 

*Many people, apparently overawed by decimals, assume without calcu- 
lation that 3.1^ must be the better approximation, presumably because 
it' is in decimal form. But notice that 22/7 - 3. 1**29, which is i 
closer to tt 3.141&) than is 3.14. 



"hard for the Greeks) to believe that we do not get all 
the points on the line (in effect, all real ilffmbers ) . 

I ■ ' 

But we don't. tt is very ea^y to construct a line 
segment uhose length is /?; for e-x^nple, the hypotenuse 
'of -an isosceles * right triangle wyth legs'of unit length.' 
I f we then place th,is segment on*the line of the ^preceding 
paragraph with its 'left end at 0, its right ^end will .be 
at /2S And /Z is irration|£, as the Greeks 'discovered to 
their consternation,* ami as you may have seen proved 
elsewhere., Therefore is not one of the points marked 
off In the construction o^f the preceding paragraph. 

^2 Decimal Representation 

Ma thematic ians now know that there are very many 
irrational numbers. In fact, as you may already\know, 
a rea^l number expressed in decimal form is rational if 
and only if it eventually becomes repeating. (A term- 
inating decimal can be regarded as a decimal that repeats 
zeros Some examples are: ' 



1 



0.25000 .'. 
0 . 5183000 
0.3333 



25 _ 1 
TW . " T 

5183 

" I070W 
i 

= T 



'o . 14285^142857 . . . = i 
- . t ' 



,0.111 



1 

5 



( v n • r i 2 

0^1333 ... * 0.333 ... - 0.2 = 4 * r = -rr - 

This is not, the same as saying that the decimal 
expansions of irrational numbers nevei; have a pattern. 
Th^ere are*some patterns which are £Ot repeating. One j # 
exampl'e is 0 ..10110011100011110000. ... . 



v *The Greek's consternation was based on much more than having their 
"common sense" jarred . Some" of their mathematical theory and even of *V 
thei p* phi losophyvds based upon the t erroneous belief that the.ratio 
of any two line segment lengths was rational. " w 
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Exercise 1 



Each oKthe following decimals can be obtained from the 

repeating decimals listed above by simple algebra. Use this 

fact to express eaph of them in the form £- , where p and q are 

i ntegers . 
a- 

(a) 
(b) 
(c) 



0.003333... 
0.766666. . . 

0. 892*857 1^28571 ^2857. . 
0.253968253968. . . 



2. THE PROOF 



* In 1767 Johann Heinrich Lamhert (1 7 2*8 - 1 777) an 
Alsatian philosopher, scientist, and mathematician, 
proved that it. is irrational (see Exercise 10). It 
follows that we- can never writ£ down a fraction or 1 
terminating decimal which equals tt ; we must always 
settle for*an approximation. ; 

v Since then ther.e have been many different proofs * 
that tt is irrational. The proof presented in this unit 
was- originally d iscovered v by Ivan Niven in 1947. It has 
the virtue that a first-year calculus student can read 
it, but it do§s -contain, unfortunately, a considerable 
. amount of computational detail. It will be much easjer 
'to follow if, before getting into that detail, we pTes^ent 
an out 1 me . <» 

> 2 . 1 ' Outlirte ^ 

^ ^ It is much easier to wo'rk computatipnally with 
the notion of rational number than with that o-f irr 
tional number. After all, ar\ irrational number is 
defined by what it is not- (not rational) rather^than 



what it is. 



So most proof that any specific number 




9 

ERLC 



1 7 



/ 



is irrational are indirect; they begin by assuming it is 
rational, and then arrive at a contradiction,. Our proof 
is no exception. / * 

If 'we assume tt is rational, then so is tt 2 , and we 
can write v 2 » p/q where 'p and q are integers. The heart 
of our proof will consist of an* intensive study of the 
express ion 

(1) K = *p n f F(x) sin *tx dx • \ 



where 

(2) _ F(x) = 



x" (1 - x) n 



i 

anywhere n is a fixec^ positive integer. Notice that K, 
depends upon ri^ (but not upon x) . • 

First, we shall make an estimate of the size. of the *' 

t . < 

integrand and, using this estimate, shall show that* for 

very large integers n, K < 1. 

1 Second, we shall perform the indicated integration, 
obtaining an expression for £ -in terms of T\0) , F(l), 
and higher derivatives of F evaluated at x = 0 an^d x = 1 , m 

*Third, vie s"hal*l derive certain properties of 

F(x) and its derivatives and apply thefo^to this exjgressio^ 

to show th^t for aj,l positive integers n , is a positive 

integer. ; ^ 

• * # 

The two italicized statements contradict each 

other. Since our <lhly assumption along the wa^ is that 

tt is rational^ this assumption must be false. 

2.2 Part One '(K < 1) * ^ - 



The function we must integrate in 'c6mputing K is 

- x) sin TT x u , ^ c ^ 

— n ~7 . ^ We are integrating from x = 0 to 

x = 1 . In this- interval ^& numerator is the product 
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of three factors, each' of which is between 0 and 1, and 
-so it jnust also be between 0* and 1. 

<T< x n (l - x) n sin Tix < 1 
0 < x " x) sin ttx < _1_ 

• 1 Jo n! n! 

We have introduced Strict inequality signs in the last 
li<ne, since .the only way we could have equality would 
► Be for the integrand to be identically 0 or identically 
^j- for 0 < x < 1. Clearly this is not so. 

. * Multiplying by 7xp n we get 

. 0 < K < =jf . . " 

This part of the proof will be complete if we can 
show that * 1 for very large*n. To show this we 

let 

. r = 2p. 

Trt^re is some positiVe integer m for which 

ee Exerc ise 2 ) , 

* • 

Then v ' 

" 1 < r! 



0 m r 
2 Tip 

and we have * P 

?p r *" „ HPS 111 + 

(r+m)! r! (r+1) (r + 2) . . . (r+m) 4 

r! r+1 r+2 r+m 

< ie! . & ; .... £ ..ml' (w n < 2 n ---L - i 

A / TT r r ^ r! z 2 m " u ' 

» y ,m times ^ ' ^ s 



Exercise 2 



Y 



Assuming p and r are given positive integers, find positive 
integer m so that 2 m > . /"Try your answer out for specific 
values of p and r 



/ naauiu i ny p ana f dre g M 

4er m so that 2 m > ^ .j'l 

" : ~ — 

* « 2.3 Part Two (Integration) ™ • 

■ ■ / 

The second part of our proof consists of performing 

the integration indicated in (1), The most direct way 

to go about this requires integration by parts and 

mathematical induction . If you have never studied 

these topics you will find aji alternate approach in 

Exercise 3 at the 'end of 'this" section. Although it is 

a little awkward it is perfectly correct. 

Let f(x) be any dif ferentsiable function, and, 
for any positive integer k,. let f/^(x) denote the k th 
derivative off(x). We shall use mathematical induction 
on n to prove that, for all non-~n.ega«tive integers n, 

(4) f 1 f ( K) sin Kx dx = miimi : (-ww 



f (2n) (1 j +- f P2n) f0) -j f, (2n+2 ) 
TZTf — + ~ TnTT . f " 



w 2n + l "Tn+I 

7T JT 



) sin ttx dx 



t 



where the upper signs applyif n is even k and the lower 
signs if n is odd. The proof, like *all .proofs involving 
mathematical induction, requires two steps: (A) Proving 



^(4) for x an initial value of n (in this case, n=0); 
(B) P/roving that if (4) holds for n=k then it holds fbY 
- . n - k+1.; - ^ 

(A) T^e first step involves 'two integrations by parts. 
Replacing n by 0 in (4) , we see tha^we must prove 
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(5) Hr(x) sin ttx dx * f C . 1 ? 1 f (°) - \L f l f"(x) sin ttx dx 

JO . 71 7T 2 JO 

v - * • • 

Jhe formula for integration by parts tell us 

f l • 1 f l 

u dv = uv - v du. ' 
Jo . 6 

"We can apply this to the left side of (5), using +* 

, *• 

u = f(x) dv = sinirxdx 

1 " 

du =' f'(x)dx . v - — cos ttx dx. 

V ' IT 



We get , . * * 

(6) [ f (x) sin ttx f - ^^cos ttx 1 + - f 1 f ' (x) cos ttx dx 

JO 7T ' 0 ^ 71 J 0 

Applying integration^ parts to the integral on 
the right, with 

■ *» *• 

i> = f ' (x) ' dv - cos ttx dx 



.v =. - sin ttx 
tt 



% du = f M (x)dx . .v =.i 

we get % - « 

[ [f • (x) cos ttx dx = f ''( x ? sin ttx * - i f 1 f M (x) sin ttx dx 

JO, 71 0 TT J o 

* » 

1 f 1 

= A - - -f"(x) sin ttx dx. 

J o _ 

i 

Substituting this*result into (6) gives us (5) and 
completes the first step in our mathematical induction. 

e t 

(B) Npw we assume (4) 'holds for n=k." In other words, 
we assume 



ii 



i r 

f(x) sin ttx 
o ' 



dx = f(i)+f«n . f"ti)*f"(o.) * 

* TT 



' rzlm 



1 



Q f v ' (x) sio tt x* dx. 

Equation (*5) , which we have just proved, applies to any 
4 dif ferentiable function f. If we aptly it to'f( 2k+2 } 
(replacing *r*(x) by ' f (2k+2) '(x) throughout)? we get ' 

V 2 * +2 ?(x)si„,xdx Ji^hn^hoi 

0 —* TT 



f J (x) am ttx dx . 

o r 



Multiplying this by + -7^77 g lv es us 
r (x) sin ttx dx 



TFT 



. 7 f (2k +' 2) (l) + f ( 2k+2 ) 



i.Zk+3 

7T 



,(2k+4), . . : 
f (x) s in tix dx . 



Making this substitution on the right side of (7) yields 



(8) 



f(x) sin tt x dx = lilllim.': f*(l)+f"(0r , 

' \ ' • 

+ f (2k) (W 2k >(o) ; 



'^ f (2k + 45 



(x) sin frx dx 



• 4 7,7"? 

whi'ch ds precisely (4) with n = k+1. Vie have comple**cf 
the proof of (4). *• . — 

Now let us appl/ this result to (1)., setting 
f(x) = F(xf). Since F(x) r is a poly'nomial'jof degree 2n, 



(xft. Sin 



12 




F (2n + 2)^ j£ identically zero, and the integral on the 
right side of (8) drojis out. Therefore * - 

(9) " P P (x) sin dx * ' 

J ° * F (2n) (mF (2n) (QK 



= F(1)+F(0) m F"(1)+F"(0) + 



' 2n + l 

7T 



Multiplying by* it and using the fact that tt 2 ^ E^- for any k, 
# (10) 'tt F(x) s in .tix dx 

• 1 ;' r(mra) 



= F(1)+F(0) 



_ F (2n) (1)+F (2n) (o y n 



Finally, multiplying (10) by p n , 

(11) K = p n (F(l)+F(0)) - p n ' 1 q(F"(l)+F"(0)) 

' ♦ ... 7 q n (F W(2n ^(l) + F (2n) (0)). • 

* This is the expression for* K promised in phase two of our 
out 1 ine ♦ 



Exercise 3 

Here is a* way to obtain (9) without using integration by parts 

V 

or mathematical «induction.. _ ^ 

/ 

K f 

Consider the function , 
g(x) = - ^F(x) cos tix + ~F' (x) b sin >ttx 

+ -V'l(x) COS TTX ** —r F 1,1 (x) Sin TTX t 

TT J TT 



-^F iv (x) cos ttx" + -^F V (x) sin ttx 

TT b F TP 



. 13 
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I 

TTSTT^M «>« ,x 1 - 5 lL 7 F< 2n+, )(x)sini l x 

where, for e^mple, F (2n) (x) denotes the (2n) th derivative of F(x) . 
(The signs change after every odd-numbered term: -++--++-,..* 
The last sign -will be + i f n is even and - if n is odd,) 

(a) , Show that g'(x) « F(x) sin ttx, 

* * * N - 

(b) Use (a) to deduce Equation (9). 

Exercise ¥ 

Compute the following integrals, 

(a) f x 2 sin Tlx dx 

(b) . 1 x 9 sin ttx dx [Do not multiply all the numbers in 

} ° i 
your answer. J 

Exercise 5 

Compute the following integrals. In each case you will have 
to make a change of variable (substitution). 

(a) j 11 y 9 sin y dy [Do not yultiply all the numbers in 



(b) , j l y sin (ny 2 ) dy 

(c) f l y 5 sin fry 2 ) dy. 
Jo f 



your answer.] 



. ' 2.4 Part Three (K is an Integer) 

What remains to be done is to show that the right 
^ side of (11) equals a positive integer. Ii: is made up 

•of sums and products of p, q, and terms of 



F (k) (0) 



and F^ '(1) for various, integers k 



the form 
> 0. Therefore 



it will be more than enough to show that Ffx) and all its 
derivatives take on integer values at x M and x = 1. 
This will show that K is an integer. That K > 0 follows^ 
immediately from (1) and (3). 

Let us first consider x r 0. We have 



already 

remarked that F(x) is really a polynomial or" degree 



14 
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tft ~ 

2n. So far any k < 2n its k derivative is a polynomial 
of 4ego-e'e 2n - k (for ^ > 2n the k th derivative is of 
course identically zero). -Further, the value of any 
polynomial at x = O'equ^ls its constant „term: - And, 
since differentiation reduces the exponent of each non- 
constant term by one, the constant term of F^(x) is / 
determined solely .by *the* "x k " tefm *of F(x).^ . •* 

v With these background facts in mind we can go to 
work. To. begin with we can expand (1 - x) n : 

(1 x) n = C 0 + C,x + C 2 x 2 + + <^ x n 

where C Q ,^ C R are integers. (It is easy to compute* 

them explicitly [see Exercise 6], but it will not be 
necessary for this proof.) . Then ' / ' 



x n (l-x n ) = C 0 x n ♦ Cx n+1 '+ Cx n+2 *'... t C x 2n 
0 » * * n 



If 0.< kj< n there is no."x k " tern in F(x), and 
therefore^ ^ (0) ^ 0, which is certainly an integer. 

-If n < k < 2n then the "x k " term is -t^ x k .' Its 
— n i 

successive derivatives are * » 

k kk-n x k-l k ( k - 1 i C k-n k-2 M^;l)(k : 2)C k , n k 

etc. Each di*f f^ent^ia.tion introduces another factor in 
front of the C k ^ By the time we differentiate* k times 
there are^k'such factors, and so we have 



ty\ k(k-l) (k-2) ..'.2-1 C, . . ' 

F (k) (0) = ^ ^n u k! c 

ni n! k-n 

k ' * (\c ^ 

Since k>n, -± is an integer, and so is F W 



^j- is an integer, and so is F^ J (0) 
(kl 

To show F v ; (l)-is an integer for each k, first 



observe that 
(12) -F(x) = t F(l - x) 



11 



r 



< It 



r 



4 



for a,Ll x. We can differentiate both sides of (12) with 
respect to x, applying the chain rule 

[f(g(x))]' -f , (g(xX)g'W 
to the right side with g(x) = L-x and f = F. We get » 

* . F*(x) = -F' (1 - x)_ 
Differentiating again (and again using the chain rule), 

F"(x) = F"(l - x) . / 
Doing this K times we see ' ^ 

F (k) (x) = + F (ki V - x) 

with the sign depending on whether k is even or odd. 
Putting x = 1* we get " 

F (k) (1) = + F (k) (0) 1 ^ 



and since we already know F^(Q) is an integer we are 
done . 



Exercise 6 



Compute the coefficients. C , ...,*C if (l-x) n 



i + c,x + c,x* + ... c;x n 

1 * A 



3. AN EXTENSION OF THE RESULT 

a As mentioned .in Section 2, the result that tt is 

iTrational is in a sense a negative result. It tells us 
what kind of ^number tt is not Another way^ of putting 
this is^jjft, if we start with the integers, no amount 
of addition, subtraction, multiplication, or division 
can produce tt exactly. But there are other operations 
we can perform on ^integers; for example, we can take 
square roots, cube rrots, etc. Is it~po$sible^ that 
one of these will produce tt exactly? . 
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o The answer is "No." This, and much more, follows 
from a result first proved by F. Lindema^n (1852-1939) 
in 1882. In these concluding section^ We shall discuss 

<-Lindemann's results briefly. ' * ^, 

v 3.1 Algebraic and Transcendental Numbers 

First we n^ed a couple of definitions. We caA a 
number algebraic if it is a root of some polynomial, 
equation with rational coef f icients^ Otherwise we-call - 
it transcendental. For example, | /is algebVaic ^smce 
it, is a root of Sx - 3 = 0, and /I! is algebraic since 
it is a root of x 2 - .2 = 0. In fact, every rational 
number is algebraic (see Exercise, 7), and every number * 
of the form % is algebraic, where m and n are positive 
integers (see -Exercise 8). Further, there are lots 'of 
polynomial equations which we haVe no idea how to solve. 
Nonetheless, we can assert that their .roots, whatever 
they may be, are algebraic. For example, the equation 

has at least One real root, and thajt root *is algebraic. 

' V/ 

: z 3 ' 1 

Exerc ise 7 

- > 

Prove that every rational number is algebraic. 

Exercise 8 

If m and n are positive integers, pro\^ that "Vn is algebraic^ 
Exercise 9 ** 

Some writers" replace the word "rational" with "integer" in 
the definition of algebraic. Prove that the two forms of the 
definition are equivalent. That is, prove both of the following: • 

(a)*if x is the root of a polynomial equation with rational 
coefficients, then it is the root of a polynomial "equation with 
integer coefficients; 
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(b) if x is the root of a polynomial equation with integer 
coefficients', then it is the root of a polynomial equation with 
rat i onal coefficients. 



3.2 lindemann 1 s Result 
V 

'''Lindemann proved in 1882 thart tt is trancendental . 
So Exercise 8 immediately tells us that, as already 
asserted, we cannot produce tt exactly by taking roots of 
integers. But it also tells us what kind of number tt is 
not (i.e., tt ^s not algebr-aic) and extends the list 
of procedures (i.e., solving polynomial equations) which 
can not give us tt exactly. Incidentally, the ancient 
Greeks could not have conceived of t^his result b.ecause 
of the strictly geometrical way in which they worked. 

ft » ; ' , 

Exercise 10 . ^ 

Lambert's original proof that tt is irrational, is based on 
the following, which he proved: if x is any, nonzero rational 
number, tJien tan x -is irrational. 

Assuming this is true, explain how it follows that if is 
i rrat i,onal . * • 
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4. HINTS FOR SELECTED EXERCISES 



3 V Use the Fundamental Theorem of Calculus. 

4. _ Use W withf(x) = (a) x 2 (b) x 9 . All derivatives of f after 

the (a) second (b) ninth are identically zero. 

' 2 2 C 

f>. (a).x = v/tt. (b) x = y . (c) x = y . Then in each case use 

(M. 

6. Use the binomial theorem. 
10. Consider x = tt . 

Jr 

5. ANSWERS TO EXE'RCISES 

. t \ J_ 1 . 1 r . v _L 2 _ 23 

" la ' 100 x 3 " 300 * D ' 10 x 3 "30 

ln-ii- 

2. Any integer greater than ^ ^ . 

* AW -i-4 ( 5 , l.M + 9^6 9-8-7-6-5-^ 91 

71 7T 3 . 71 7T 3 1 7T 5 ^* 7T 7 7T 9 

5 . (a) ,.0(1.9^^9^.9^7^,9!, 

& * it* 7T 5 7T 7 7T 9 

W ? • " ' 

i \ 1 2 



6. C, = (-1)' 



n] _ (H) 1 n! 
■ i! (n-i)! 



7. r is a solution of x - - = 0. 

, q q 

8. r Vn is a solution of xf 1 - n - 0. 

9. (a) Given a polynomial equation with rational coefficients, 

multiply it by the least common multiple of the 
denominators. The new equation has integer coefficients 
and the same roots. , 
(b) Any polynomial wi th* integer, coeff ic ients is already 
a polynomial with rational coefficients. * 

10. Take x ■ tt. ^ I f tt were^ rat i onal , then tan x - 0 would be 
i rrationa! . 
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1. HISTORICAL BACKGROUND . 

t - . 

1 . 1 introduction , t ' » 

The number 77 is defined as the ratio of .the 

circumference to the diameter of a circle. It is 

of central importance in innumerable mathematical 

and scieAt>fic results,, including many which do not- 

appear to have .anything to dcf with circles. 
> . ' X 

For thousands years^jnathematicians and others 
have been fascinated by the challenge of- determining 
the value of^r as precisely'as possible, and now 
computers have gotten into the act. In T967 a com- 
puter *in Paris determined 7i t to 500,000 decimal places. 

It is undeniably true that computers, because - 
of the phenomenal speed with wljich they can do arith- 
metic, have" enabled us to obtain "many more digits t 
of 71 than ever before. In 194& (j-ust before computers) 

v the record was 808 v places. But all „the machines can 
do is use formulas supplie&flto them by us slowpoke 
humans , ^who^re much more intelligent n% matter 
what-you may ha've heard 1 tp the contrary. It has been 
said of the computer that "the most intelligent thing 

• it is capable of doing withou't the help of its pro- 
grammers is y to go on strike when required *to* work 
without air conditioning,"** ^ 

Well, then, where' do .these formula's 4me from? 
How do "mere" humans, find the /value of 77, even without 
the aid oi computers? Afte? all, 808 places isn't 
peanuts. 'Certainly no one could measure the circum- 
ference of any circle, however large ,< with- that 
accuracx. 

Beckmann, page 102. 

* ' ' \ 
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' 1.2 Early Circle Measurers 

* About 4000 years ago the Babylonians, who knew 
nothing about It except its definition, probably did 
actually measure circles. They may have, for example, 
drawn a large circle on tfie ground and marked off its 
diameter on a~rope. They could have "then observed 

.^that laying tfris rope segment along the circle three 
times almost^ but not quite, tpok them around the 
circle.* Perhaps they then measured the excess with 
another rope and determined that this shorter segment 
went mto^the original segment 4 about eight times (we 
- now know seven would have been more accurate). At 
any rate, they came up with the estimate tt * 3^ . 

1^3 Archimedes and other Polygon Mea,sWers 

* In about 240 BC "the Greek mathematician ArchrTnedes 
(287-212 BC) became involved in our/^ale.. Archimedes 
is universally regarded as one of the two or three 
greatest men in the history of^estVrn mathematics. 

I He was apparently the first person who attempted 
to estimate Ti.in any wa^ other than literally 
measuring the c ircumf erence- of a circle. Instead 
he said,\in effect: "Let mfcstart-with a circle of 
^diameter erne'. Its circumf erence will of "course be 
tt . Instead olF trying.' t;o measure ^the circumference, 
I'll" inscribe a polygon and. c ireum'ser ibe another 
polygon. Clearly the circumference tt is between 
the two perimeters: If I choose the polygonswisely , 
I should be able to use the geometry I know to compute 
these per imeters . " *" 9 

- For an illustration of Archimedes' thoughts look 
at Figure 1, whi.ch shows a circle of diameter one 
with inscribed and clorcumscr ji)ed regular hexagons., 
Each of the hexagons consists of six equilateral 
triangles, so the calculations are straightforward. s 
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Figure 1. Archimedes' calculation of v. 
' *Nie circle has diameter one. 

For^ the inscribed hexagon AB = OA = y, so the perimeter 
u - 6 *,a£ = 6 x i ■ 3.* For the circumscribed hexagon 
we can' observe that 0F £ = > so 

2 _ 1 ' " 



5F = j 



But also, 



0F®= mi 2 - UF 2 = TJU 2 - [\ CE) 2 ' 
'/ = Uc 2 - (± uT) 2 ' = J JC 2 . ' 



1 



Equating these two expressions for OF 2 we ge£ 

3 777*2 _ 1 



so 



. and 
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Then the perimeter 

P = 6 x CD = 6 x OC = 2/3 . 

We could have used trigonometric functions, since, 
we know the*value of sin 60°, but they were introduced 
long after Archimed^',, and we wanted to show that he • 
did not really need them, , 

We've shown 

3 < 7T < 2/3 • 

.or . " 

0 

3.00 < TT < 3.47; - 

This result in itself would not be worth the 
trouble, but Archimedes used some rather clever 
geometry to compute the perime'tors of inscribed 
and circumscribed regular 'polygdns with 96 sides, 
,He got - m 

<•■<>{ 

or, 'in modern decimal notation, 

\ - 3.140 < tt < v 3.143. 

As a modern mathematician would put it, he determined 
tt to two decimal" places. 

Although Archimedes brojce away from di«rect 

measurement of circumferafces, his method of estimating 

ti was still driectly related to the fact that tt is ^ 

the Circumference of a circle with diameter one. jsJ . 
C th 

Until the 17 century all attempts to compute ti 
were based upon this fact, 6r upon the closely 
related fact that such a circle has area tt/4. As * 
of 1630 the record for digits of ti was apparently 
35, set by an otherwise obscure Dutch mathematician 
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named Ludolph von Ceulen (1539 - 1610) who is supposed 
to have based his work on polygons having 2 62 
.(? 4.6 x*i0 18 ) sides and to have devcfted most of 
his life to this task, 

1.4 Analytic Attempts 

t h * 

. In the 17 century the search for digits of tt 
began to take a fundamentally different tack for the 
first time s^nce Archimedes. Mathematicians began 
to turn away from circles, polygons, and other geometric 
considerations. The mathematical concept which were 
to become» M Calculus" under Isaac Newton (1642 - 1727) 
and Gottfried Wilhelm Leibniz (1^46 -1716) were in 
their infancy then, -and mathematicians were just 
beginning to understand the notion of "sequence of 
rational numbers.'" From the 17 th century on, all 
attempts to approximate n , with iyr without a computer, 
have amounted to finding sequences o£. rational numbers 
whose limit is tt. It was* this approach which made it 
possible to extend the rejpord from 35 digits in- 1630 
to 527 digits in 18f4 (before desk calculators) and 
808 digits in 1948 (before computers) , 

In this unit we shall derive one of the earliest ~ 
of these sequences. It was originally discovered by 
the English mathematician Jo£n Wallis (1616 - 1703) 
inabout 1650. 

<*> 

2. THE ■ WALLIS FORMULA 

2.1 Outline; D efinition of I 
' ji 

Our derivation of the Wallis formula is, based upon 
a study of the definite integral 

sm n x dx , 

0 

5 

i 

9 o 
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where n can be any positive integer, or zero. 1 Since 
we will be talking about thi| integral quite ajlot, 
let's give it a name: 



■■•J 



tt/2 



sin x dx 



We'll proceed in five steps. 

1) We'll use integration by. parts to obtain 
a formula expressing I R in terms of I n-2 



2) 
3) 

4) 



for n ^> 2. 
We'll *£c this formula to compute I 



n 



We'll then compute I n+1 /I . It will turn 
but to be - times a certain rational number 
r^ which depends on n. . 

We'll go back to the definition of i' and 
sfew direptly from it, without using (1) , 1 
(2) or (3) , that ' 



1 im 



n+1 



T 



= 1, 



2.2 



5) Combining (3) and (4), we'll observe that 
2 

_ x r n -^1 and therefore 2r n + n , 

{2r } = {2r n , 2r ? , 2r_, 
n 1 2 3 

sequence . 
Reduction Formula for I„ 



The sequence 
. . . } is the Wallis 



We'll begin by attempting to compute I using 
integration by parts. Using the formula. 

ju dv = uv - jv du 



with 



u * sin^x 



du 



(n-1) sin n j|j x cos^x/ix 



! 

dv = sin x dx 
v = -cos x 
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we get 



sin n xdx = -sin"" 1 x cos x + (n-1) |sin n " 2 x cos 2 x dx. 

Replacing cos 2 x by 1 - sin 2 x" in the last term: 

(1)* jsin n x dx •* 

* n-1 
' = -sin t # cos x + (n-1) 

' - (n-1) Jsin n x dx . 4 



. .n-2 , 
sin x dx 



Novtf evaluating from x = 0 to x = | : 

ir/2 

*(n-l)I n _ 2 - (n-1)!^ 



t . n-1 

I = -sin a x cos x 

n • , 5 



0 



I n = (n-1) I n _ 2 - (n-1) I n 



(2) 



« . n n- 2 . 



J 



Exer-cise 1 . From (1) obtain an iteration formula for /sin xdx, and 
use it to determine: * 



(a) /sin 5 x dx; (b) /sin 1 * x dx. 



*g . 3 Computation of I R 

If you have never seen a derivation like the one 

leading to Equation (2) youmay think we have failed 

*in our attempt to compute I using integration by 

parts. After all, we have "nierely" expressed one 

unknown 'integral in terms of another. * But I 

involves *a lower power of sinx than does L . Formula 
r n 

(2) is called a reduction formula. If we start with 



the direct computation 



Jo 



sin 



rTi/2 

x dx = 1 dx'= x 

1CI * i 



Till 

0 



TT 
1 



we can use (2) again and again to work our way up to 
I f&r any even n: t * . i 
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n - 2: 
n = 4 : 
n = 6: 



1 

1 • 


»0 


_ 1 
" I 


TT 

7 


3 

T * 


! 2 


3 
T 


1 
7 


5 


! 4 


5 

" 5" ' 


3 
T 



TT 

7 
l 

"7 



TT 

7 



and in general 
(3) I 



'n-1 
n 



n-3 
n-2 



for any even n. (This is called an iterative* process, 
from the verb iterate to repeat). 



If we 


start 


with 






















7T/2 t 






h 


-\ 

0 


sin x dx 


1 


COS X 


= 0 - 

0 


(-1) 


obtain: 






V 










n = 3: 


' h 


2 
3 




1 








n = 5: 




4 

5 




2 
3 


• 1 






n = 7: 


h 


_ 6 

7 




4 

5 


. 2 
3 


1 





and in general 
for any odd n . 



n-1 . n-3 
n n-2 



2 S 4 



Computation of I n+ f^I 



Now we are "ready to consider the ratio I n+1 /I . 
The formula for this ratio depends 'on whether n is even 
or odd. Let us suppose n is even. (It will turn' out 
that this is the only case we f ll have to consider.) 
Then (3) gives us the demoninator l^. Since n is even, 
n+1 is odd, and we can obtain, the numerator I n+J from (4) 
if we replace" n by n + 1 in that formula: 
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r c> t a n m n-2 6 4 2 - 

- (5) ^n+T rTT 4 4 4 7 ' 5 ' T 

| £or any eyen n. 

Dividing (5) by (3) : 

I n n-2 6 '4 2 n 

n+1 _ n+T n-1 ' ' ' 7 * I ' T ' 1 



< 



T" n^l n^5 5 3 T 



n 



n ' TPT ! ' ' !>•* T * I * 1 



Taking the factors alternately from the "numerator and 
denominator, and remembering to invert those from the 
denominator , ; 

*n + l _ _n_ _n_ t nj_2 . n-_2 6 6 4 4 2 2 

I n ~ n+1 ' n-1 ' n-1 n-3 ' ' * 7*5*5' 3 ' 3 ' T 

for any e^ven n. It is conventional to write these factors 
in the opposite order: e "* 

4 

(6) *n + l .2,224466 n n 

~T" ? 1 T'T'S'T'T'T"' n 7 ! * n+T 4 



2 . 5 The Limit of I 

n+ 1 n 

Let's calculate the right side of (6) for the 

first'few even values of n* 
* 



'2 



» - 4- -£ 1.^2.2.44 a 128 0. 9054 

1 I^-tt 1,3 3 5 • 4 5tt 

_ * c. a *7 w 2 , 2 . 2,4 4 6 6 4608 ^ n n , 17 

n - 6. T-.y.l^.y.y.y.y^ - . 0.9313 



6 

Io 



0.9461 



n - 10: yii = ... * z 0.9556' 



L 10 
! 13 

n = 12r = ... z 0.9623 
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It looks like this sequence may be approaching one. 
Can we prove this? 

We can prove even more. As^stated \n Section 
2.1, ue can go right back to the definition of I • 
and sho> thzyt ' 

I hn *n+l . 
« n**> I 

n 
«<* 

In other words, the sequence ^^l' 1 *> ' l 4 /l ~> ■•■ 
approaches one. Since our sequence L/L, U/I,, 
U/I^, ... consists of every other term in this 
sequence, it must also approach one. 

Remember that all of the I are integrals from 
x = 0 to x = For 0 < x < ~ ue know that 

0 < sin x < I . Since 0 2 = 0, l 2 * \ , a 'nd y 2 < 
all y between 0 and I ue ha\ e 

(") 0 < s in 2 x < s in x <_ 1 . 

Since sin x > 0, sin n x > 0 for all integei s n. Ue 
can multiply (7) b> Mn n x, getting » 

0 < *in n+2 x < sin nM -x < sin n x 

and therefore * 

- 0 < .1 , < I < i . 
- \n+2 - n+1 - n 

(If 0 < f(x) < g(x) throughout an interval [a>, b], 
, then 

b b v 

f(x) dx < I g( X ) dx , 
a , a 

since these in tegrals represent are^s, one of which 
is contained in the other. See Figure 2.) 
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f (X) dx 



j b g(x) dx 



Figure 2. Graphical proof that Tf (x) dx < f b g(x) dx when 0< f(x) ^ g(x). 

J a ~ ' a ~ " 

Since I R * > 0 we can divide by it*. 



I a0 I . 
0 < J^l < "ii < 1 



- I 



n n 

or, using (2), with n replaced by n+2;^ 

n+1 - *n+l , . 

n+7 - "1 - * 
n ■ 



Now suppose rw«>. ^WreTT certainly -> 1. 



The middle expression Is trapped/6etween two expressions 
which are near one whenW---T^Targe (one of them 
actually; equals one) 1 , /Ad sq. it too is near* one when 
n is large. That is(j\t approaches one as n approaches 
infinity : 



1 im n+1 



1 . 



(This is sometimes* called the squeeze principle. See 
Figure 3.)v> * 

v 11' 
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Figure 3. The squeeze principle for {£±1} , U^/I,,}, and 
= {1,1,1,...}. 



2.6 The Wallis Formula for' tt 
— - . 

,We have now proved what we guessed might be true 

earlj^in Section 2.5; that the sequence whose first 

six terms we calculated there does indeed approach 

one. Therefore its terms must eventually be near. one. 

A bit more precisely, 



(8) 



2 2 2 *4 4 6 6 
tt * T* 3 "I' S " 5 "7 



? 1, 



n-1 * fFT 

-v. 

and the left side can be made as near to one as we like 
by taking n large enough' (and even). 

Multiplying (8) by we get 
fy large N even n This is the Wallis formula. 
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3. "exercises 



(a) Show that 2-4 «6i • • (2n) = 2 n n! for any positive 
- - - -^ere n! denotes l-2«3"-n. 



i n t eg er n , w^c 



(b) Show that 

2 -2-4 -4- ■■(2n)«(2n) 
1 . 5. -s- (2n-l) • (2n + lT 



(n ! ) 2 



4n 



[(2n) !] 2 (2n+l) 



for any positive integer n. 
(c) Show tha\ 

lim (n!) 2 2 2n 



(2n)! 



This problem presents a derivation of>another 
formula for it. The formula is credited to Leibniz. 

(a) Show that 

2n+2 



J - x 2 + x k - X 6 + . . . + x 



2n 



1 



1+x 2 1+x 2 

for any real number x a*hd any even positive 
integer n. 

(b) Prom (a) deduce that 

ill > i f 1 



1+X 2 



Cc) Show that 

r 



Y 



■ 0 < f 1 *^dx < f^m*? dx u l 
• " J 1+x 2 " i 2n+ ~ 



0 A " ~ 0 
(d) Finally, deduce that 



lim 4 h 1 + 1 1 .... + 
n+<» j 5/ 7 2n+3 



This is the l.eibni'z formula for/ tt . 
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4\ (a) Integrate the formula m 3(a) from 0 to u ' 
J* (where 0 < u < 1) and then, by steps similar 

♦ to tho^e in problem 3, show that ♦ * 

where n is a large even posi'-tlve integer. 

(b) What value of u should give ^ on the left < 
side of f(9)? * 

(jc) Use this value of u.to obtain ara5T?\er 
formula for tt . - ► 4 




The Project would like to thank Charles Votaw oi Fort 
Hays State University and Solomon Garfunkel of the University 
of Connecticut for the-ir reviews-, and all others who assisted 
in the production of* this unit. " * 

This unit was field-tested and/or student reviewed in 
preliminary form at Montgomery College, Rockville, Maryland; 
^ \}suMeT9ity of Louisville, Louisville. Kentucky, Lycoming 
- V College, Wi 1 1 iamsport/Pennsylvania; Fprt Lewis College. 

Ourango, Colorado.;, afid Franklin College, Franklin, Indiana, . 
and has been revis/d on the basis of data received from 
these sues. * * « 



0 

36 



4. HINTS AND SOLU TIONS TO EXERCISES 

' ! 



(aj j sin" [x dx = ■- SinT1 " 1 n X C ° S X ♦ ^ j sin n_2 
Starting with j sin x tlx s • cos \- 

n = % ■ LnJvHv - - srrrx cos x 2 
« n - ^ . j sin x ox - — * - j cos x 



- S: fsin'xdx = :_si "X cos ' x + 4 f s , n , x ax 



1 ..»_... . 4 



1 



= "jSin x cos x - y$ sjn 2 x cos x 



COS X. 
1 ^ \ 



(b) n -I 2: Jsin 2 x dx = -sinx^osx + \ jl dx' 



2" s in *3*sos x + j 



n = 4: (sin^xdx - iiJl!*.^* + 5 j sin 2 ?dx 



1 3 

- js in x, cos x 

3 . 3 
- o- s in x ccs x + ~x . 



(a) 2-4-6---(2n) = (2 - 1 ) ( 2 . 2) (2 - 3) . . . ( 2 ' nj 
• <» 

= f 2;2;2---2 )(l-2-3---n) = 2%!. 
„ ' i ' 

n 

(b) " Multiply "loft side by jj'l'l" 44^444 . 

» * 2 • 2 • 4 - ^ . . . (2n) (2nJ 

Show jthat the* new numerat or is (2 • 4 • 6 - « . ( 2n) ) k 
and use part_ (a) . 
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(c) The Wallis formula together with part (b) 

says 

„ , 2 (n')^2 4 " 



[(2n)!] i (2n+l) 



Taking square roots we get^/ 
(2n)! /2n+l) 



or 



i lllB (n!) 2 2 2r V? / - 
f 1 1 m -- — m — ■ = /tt • 

n-°° (2n) ! »'7n+T 
Multiply by /n//n and observe that 

/In 



/In = /~7TT 
/2n+T / 2n+1 



■* 1 as n-*°° . 



(a) Use the fact that , 
x 2n+2 = (l + .x 2 j(l-x 2 ^ 4 -x 6 + .V. + k x 2n ), 

which can be confirmed by dinrct multiplication. 

(b) integrate the-equation in 3(a) from x=0 to x = l: 



f 



2n+l 



'l 



tan 1 x 



1 x 2n + 2 



dx . 



0 0 



1+x 2 



2 2 
(c) Since x > 0, l+x» > 1, and therefore 



(x 2 " +2 K 
\ =r < x 



2n + 2 



Recall that if f(x) l'g(x) for 



all x between (T and 1, then 



1 



f(x)dx < 



1 

g(x)dx. 



(dj Yn+rj "* 0 as n-**. Therefore, using the result in 
3(c) and the squeeze principle, 
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f 2n + 2 
7 



dx - 0. 



1 + x 



Applying this to the equation in 3(b) we get 



i 1 ♦ 1 1 ♦ 

1 "5 T> 7 



Taking u = — : 
/3 



1 1 



1 



1 



/I 3-3/3 5-3 2 /3 7-3 3 /3 



^ (2n+l)3 n /3 



_6_ 
✓3 

2/3 



1 - * ♦ 1 

j ■ 3 



1 . 1 



1 



5-3 2 7-3 3 
1 



5-3 2 -7-3 2 



1 



(2n+l)3 
1 



(2n+l)3' 




\ 
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1. INTRODUCTION * < 

1.1 The Experiment , 

Certainly you have met the number tt'(z 3.141592) 
very frequently in your study of mathematics. It is 
defined as the *ratio of the circumference of a circle 
to its diameter, but it arises in many places which 
appear to have nothing whatever to do with* circles. In « 
this unit we shall describe an experiment where n 
<does indeed arise unexpectedly. \ I»t provides a "fun and 
games-" method of approximating tt . 

You can easily set\ip and perform this experiment* 
'yourself. /All you need is a large flat tabletop, enough 
paper td cover ^t, and some thin object a few inches 
long, such as a toothpick. As, the unit title suggests, 
Comte de Buffon (French; U07-1788) referred to a needle 
when he first discussed this experiment. Out of respect 
for the Count, we shall refer t& the object as a needle 
throughout this' unit. A toothpick, however, is probably 
much more convenient and somewhat less dangerous. 

Cover the tabletop with paper. On this paper draw * 
a bunch of parallel lines. Make sure the distance from 
each line to the next is exactly equal to the length of 
the needle. 

Now pick up the needle and toss it onto the table. 
When it comes tb rest it will either cross one of the 
lines you have drawn, or it will not. Let us call this 
"toss* a "success** if the needle crosses a line. ~TTclT 
• up the^, needle and toss it again. Continue doing thisV 
keeping track of the total numbers of successes and 
the total, number T of tosses. 
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1.2 A Sample Run 



Let us consider the ratio S/T; that is, the fraction 
of tosses which are successes. A 'running calculation of 



S/T 


might look like 


this 


» # 


- 


\ 




* 






* 

TABLE I 






T " 


S 


S/T 




4 

T 


c 

0 


r It 

S/T 


1 


1 


1.00000 




e 

102 


65 


0.63725 


2 


1 


0.50000 




103 


'66 


0.64078 


3 


1 t 


0.33333 




10*t 


66 


0.63^62 




2 


0.50000 




105 


67 


0.63810 


5 


3 


0.60000 










6 


^ 


0.66667 










7 


5 


0.571^3 






6368 


0.63712 


8 


5 


0.62500 




9996 


6369 


0.657*5 










9997 


6370 


0.63719 










9993 


6370 


a.63713 


100 


6*t 


0.6^000 




9999 


6370 


0.63706 


101 

r 


< 6*t / 


0.63366 




10000 


6371 


0.63710 


In practice 


it would 


be 


unreasonable to 


toss the 


needle 


10„ 


000 times 


(although in Section 6 we 


discuss how a 



computer could easily • simulate *this , and much more). 



But we want to make a point. Namely, the values of 
S/T at first fluctuate wildly, but after a very large 
number of tosses they tend tb settle down. We shall use 

calculus to show that in a certain sense they are most 

2 

likely tt) settle down to a value near — . To understand 
exactly what^hat means, w,e must make a brief digression 
into probability. , 



Exercise 1 

Perform the experiment described. Toss the needle i 00 times. 

S 2 * 

Then write z — , where S is the number of successes you have 

recorded, and compute an estimate of tt. Finally, compute the 

percentage error in this estimate. 
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2. PROBABILITY 

2.1 Equally Likely Events 

Lf we toss a coin, it may come down heads or tails. 

Assuming the coin is honest, there is every reas'oh in 

the world to think these two possibilities a/e equally 

liie ly. We* say' that the probability of heads is j and 

1 

the probability of tails is 

If we throw a die, it may come up 1, 2, 3, 4, 5, 
*or 6* Again, these possibilities ai^ equally likely. 
We say that the probability of each o£\them is ^. 

Generally, if there are t possible results-, of which 

exactly one win. happen, and if these results are equally 

likely* we say tn>t^feach of the , results has probability 
1_ 
t* 



Exercise 2 

If you pick one card from a welj--shu^fled deck of playing 
cards, frnd the probability it will be the Jackof Hearts ? 



Exercise 3 

Suppose you toss a coin with your left hand while throwing 
a die with your right hand, and then record the combined result; 
for example, Heads - 6 or Tails - 3. Find the probability you 
will *get Heads on the coin and k on the die. 

Exercise A 

Suppose you throw one die with y 0ur left hand while throwing 
another die with your right hand, and then record the combined 
result. For example, you may get Left - 3> Right - 5, which you 
could abbreviate (3,5). *0r Left'- 5, Right -3, abbreviated (5,3). 
Find: , * , 
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j y 

(a) the probability you will get k on the left die and 6 
on the right die, 

(b) the probability you will get "boxcars" (6 on both 
dice) . 

\ —. 

2.2 "In The Long" Run" ^ 

There is another way of interpreting probability. 
Getting back to the die as an example, suppose we were 
to throw it 10,000 times. Since the six possibilities, 
are equally likely, we would expect that each of them 
would come up about the s£me number of times about 
^ of the 10 ,000 tosses woul<give 1, about g- of them 
would give 2, etc. We' can .say that the probability ^ 
is then a prediction of about what fraction of t x he tosses 
will give a certain result. 

Sometimes probability is defined in terms of this 
kind of prediction of what will happen "in the long run", 
rather than how likely an event* <j.s* on a "one-shot" basis. 
We must be careful, though. This "lon^ run" view o*f 
probability is not, and-cannot be^an assertion of 
exactly what will happen. For one^thing, ^ of 10,000 
is 1666j, which is not an integerJ So*we certainly 
cannot get exactly that number of* y say, 5 f s. But it 
would also be wrong to interpret the probability as an 
assertion that exactly 1666 or 1667 tpsses will result 
in a 5. The die has no memory* and can't keep count. 
What it does on each toss is not determined by the pre^ 
vious tosses. If it comers up S^on one toss* it> 'cannot , 
say, "Well, I|d better lay off »5 for the next* six tosses^ 
If it fails tb come up 5 for five or six consecutive 
tosses, it cannot say f \ f IHey I f m overdue. Better make 
it a 5 this time." And so it^raay'come up 5 a bit more * 
or less than predicted. It may even run a string of t£n 
consecutive 5 f s, although this is' extremely , unlikely . 




To impress this* upon your memory, think again of 
tossing a'cc^tn. The probability of getting heads is y. 
In the long run, about j of all tosses will be' heads. 
But this does not mean 'that if we tos»s a coin twice, 
exactly one of the tosses must result in heads. The 
.most likely number of heads is on^, but we would not 
be the least bit surprised *ttf^ see it land heads both 
times. If it happened to land heads the first time, 
we would certairily~n\>t say it was guaranteed tc5\ land 
tails the next ^ime. 

2.5 Compound Events 

Now suppose we toss a die once, and are interested 
^ how likely it is to come Up greater than J (that is, 

5 or 6). It 'Should come up 5 about of thWtime, and 

6 about ^ of 'the^Bime . This adds up t^^/or j, of the 
time it comes u^Pgreater i^han 4. We say the probability 
of this occurring is ^ (or i) . 

Generally, if there are t possible equally likely 
results, of which exactly one will, occur, and if s of 
these satisfy a certain Condition, the -probability that 
this condition will be satisfied is 

. r 

Exercise 5 , » » 

If yoi/throw one die, find the probability it will be even. 
, Exercf se'6 

.In the experiment of Exercise 2, find the probability of 
each of the following events:' . 

(a) ttie card Is a Jack 

(b) v ,the r card Is a Heart 

(c) the card is /a picture card 

(d) the c£rd Is an even-numbered card. 
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Exercise 7 



(c) 



3. 



the sumbf the two numbers will equal ten, 
the number on the right die will exceed the number 



In the experiment of Exercise 4, find the probability that: 
(a) 

(by 

'on th£'le?t die 

the tjto numbers will have an odd sum and* an odd 
product. 



THE THEORETICAL RESULT OF THE EXPERIMENT 



3.1 Statement of The Resul t * 

1 " % 

Now we dan get back to that needle' on the tabletop 
and state precisely what we are going to prove, and how 
it is helpful in approximating v. We'll prove that the 
.probability ihe 'needle, will cross a line is K That is, 
in a very lafge number of tosses, about | (50.63662) of' 
them wil'l be "successes." The r.esult i, independent of 
the length off "the needle. ■ 

hi Application of the Result 

count T the aPP i y K thiS reSUU ' t0SS ne6dle « times -d 
count the nUr s of successes, the faction § should 

be close to/-, although our discussion in Section 2 2 , 

suggests thkt you shouldn't be too optimistic 

degree of Accuracy. Then write'*- . I and „ solve ' ^ 

In the* example of Table I , t ■ 10,0,00, s = 6371 
TJ:M S7 1: J"' 5 ^ 16SS than °- 1? *« th. true' 

■ ~ 3 r • " d S ° 1Vlng for * 8 ives 

- s.iivm (again within jQ. in 
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Exercise 8 * 

Show^that with t » 10,000 it is impossible to get better 
than three decimal place accuracy in, estimating tt with this 
experiment. 



4. PROOF OF THE RESULT 

* 

4.1 Locating the Needle Numerically 

Although the application of our result depends upon 
the "long run' 1 interpretation of probability, it is more 
convenient to' use the "equally likely events" interpreta- 
tion in proving ^the result. The events will be the various 
positions in which the needle may land. To base any cal- 
culations upon the needle's landing place we must first 
decide upon a scheme for describing the landing place 
numerically. , 




Let us assume """the lines on tfye paper run east-west^ 
and that the distance between consecutive lines is L. 
*A little thought sho'uld^convince you that two numbers 
will tell >^)u all you need to know to decide whether or 
not the needle is on a line (see Figure 1) . The] first 
of these (call it y) is the distance of the southern 
end of the needle from the nearest line to the south. 
(If the needle's southern end should happen to be on a 
line, set y = 0. If the needle should happen to lie 
east-west, think' of its western end as the "southern" 
end.) The second number (call it 8 is the angle the 
needle makes with a ray running eastward from its southern 
end. So we have an ordered pair (y,8) with 0 < y < L 
and 0 < 6 < tt v 

Notice that we have simplified matters by restricting 
our attention narrowly to what concerns us- The ordered 
*p#ir does not really tell us where^ the needle is (how 
far east^ii it? which line is it straddling?) but, as / 
we h*ave S'aid, it does tell us whether the 'needle is on 
a line. In fact, you can see from Figure^l that this 
happens if and only if 

(1) ( y + L ,sin 6 > L. 

* » 
4.2 Equally Likely Positions :"- 

A problem arises when^we attempt to list the "equally 
likely events" the possible positions' of the needle. 
There is no difficulty finding them; the trouble is that 
tliere are too many ^of them. The southern end of the 
needle is just as^ likely to be^anywhere as anywhere else. 
The n%edle is just as'likely to be oriented in any direc- 
tion as in any other. In other words > all possible pairs 
\y,9), with 0\< y < L and 0 £ 0 < tt, are equally likely. 
'But there are infinitely many of them, and an •infinite 
subset of these satisfy (1^., Our definition in Section 2.3 



51 



8 



<• applies only to finite situations; We can make nothing 

* ~m + ~ . CD ~ 

of the ratio 



4*3 A Finite Approximation 

Let us replace our problem by one which involves a 
large but finite number of equally likely events, and 
wjiose answer will be a good approximation to that of 
the original problem. Let us pick a large positive 
integer n and partition the interval 0 < y < L into 
n intervals: ~~ 

M 

(2) y = i v ' = — v = iii v _nL_ T 

Similarly, we sh^ll partition the interval 0 < 6 < tt 
into n equal intervals: 

(j) \-l e 2 .2jr, .... e^if e n = S2L=*. 



0 *(n-l)n 

1 7r/n 2n/n 3n/n n ; 


( 


> i i 

l X (n-l)L , 
' L/n 2L/n n 1 

1 1 1 1 


1 i 



lere *are now n 2 pair! 6f the form (y i ,6j) wftfere i and 
j a^e each integers from 1 to n inclusive (it may happen 
that i* = j ) . ' 

Now let us imagine that when we throw the needle, 
instead of recording its actual position (y,6), we 
record (y i> Q-) 9 where y. and 6^. are the smallest numbers 
i n (2) - and f3) greater than or equal to y and 9 respec- 
tively. This amounts to pretending the needle is slightly 
north of its actual position, and rotated slightly 
"counterclockwise . 

In Figure 2 we have illustrated this for the case 
n » 8. The needle has fallen in the position (y,6) where 
-g- < y<-j- and -j- < 6 < ^ . Therefore we replace y 

AT 7 ' » — 

by y t ; T and e by 8, = f , recording (y,,,e 3 ) as the 
approximate position of the needle. 1 

ERIC $2 



« 



Approximate pos 



True position of needle (y,9). 




Figure 2. Typical approximation-of ^needle position (n « 8) . 



We shall define p R to be the probability th<5 needle 
crosses a line after it has "been moved in this way. This 
amounts to the probability that 

(4) y. + L sin 6. > m L. 

l » ' ^ 

Since, for any actual position (y,8) of the needle, 
Yi 2 y and 8^ s 6, and since these approximations become 
better and better as n + », it follows that p R approaches 
the true probability which we seek, The result stated in 
Section 3.1 can thus b,e expressed 



lim p « £ 
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4,4 Counting Successes in the Finite Case 

' - Since the intervals 0 < y £ L and 0 < 6 < n are 

.each partitioned into equal subintervals , e.ach of the n 2 t 

ordered pairs (y^G.) is equally likely. Therefore, ' i 

according to Section 2.3, all th'at is necessary to 

determine p is to count the number of ordered pairs ** 
y V ' m 

which satisfy (4). If this numher is m , then p = — . 

n n 2 , 

.* " n . 

To determine iri n we begin by, looking aj^any one 

particular S.TjrThere are n paris involving this 6., Of* 

J« * 3 » 

these, the ones satisfying (4) are those for which * 



y. > l (r- sin e^) 



iL 



n > L (1 >sin 



(5) 



^ i .> n (1 - sin i^) 



The right side of (5) is between 0 and n. Let k. be the 

% * * * 3 

integer part of the right side; that is, the unique integer 

such that 

(6) - k.. < n(l - sin < k.. + # 1. 

' If kj = 0 then n(l - sirt is between 0 and 1 .and 
therefore less than any positive integer . 0 S»- ( 5) is 
true for all permissible values &f i (i = 1, 2, . .., n)'. 
If kj f 0 then it is a positive integer, and (5)'is% 
falee for all values of i up .to k^ (i f 1, 2, . ..,Jkj) 
and true for all values of i aftfr that (i = k. + l ~ 
kj + 2,* 1 ,.,, n) . In either case kj counts the Inijmbei' of 
ordered* pairs' with this particular §. for whiGh (5) is 

Jalee, and theTefore n - kj counts the number for which 
it /is true. .So * 



(7) 



% - I (n -'k,) = n 2 . - I k. . 
j-1 - 3 ~ j-1 3 ' 



Figure 3 and Table 2 illustrate the determination 
of kj and of m n for n -'7. The figure shows clearly 
th*fc, for each j, kj denotes the largest' i for which 
the needle fails to-cro$s a line. * 

; 1 ■ 

Exercise 3 

Find all the k- 's^tj r 1 . 2, 3,* ■ ■ • , ti) , m , and- p : 
J s. r n n 

(a) for n = 10 * 

(b) for n = 100 (this is easiest if you write a -little 

computer program) t r ' «» 

Some of your computations ^i \\ be helpful in Exercise 10. 

— — «— 

4.5 Approximating the Count 

The result in (7) easily provides the value of p 
for any .particular n, but not in a form convenient for 
computing the limit. To achieve that end, we begin by 
noting that (6) tells us^ 

k * k. = n (1 - sin li) - 

with an error less than-one. This approximation leads 
* * 1 1 

to an error of less than — in the value of p , as you 

n *n ' J 

should be able to* verify (see Exercise 11). It follows 
that 

t 

n k. z rt - n(l - si* 

3 f n ' 

' ri - k. = n sin'il 
• 3 . n 

and 'this is approximately the number of ordered pairs, 
with this particular 8^ , for which (,5) .is true. 

If we make 'the same approximation for each 6j and 
then add the results, we obtain 

CO 

m n sin , 

n y?i n . 

12 

' ' . .55 



. TABLE II 

Calculation of m n (Section 4,4) and of the 1 
approximation of ^Section 4,5) for n=7 



j 


n 


ill 

sin J — 
n 


n (1 - sirr 

n 


k. 
J 


n - k >. ' 
- J 


n &i r> 

n 


1 


IT 

7 


0.4339 


* 3.9623 


3 




3.0372 


2 


2tt 
7 


O.78I8 


1.5272 


9 1 


' *6 


5.1*728 


3 


3n 

T 


0.97*9 


0.1755 


■ 0 


7 


^6.821*5 


k 


7 


0.97^*9 


* 04755 . 
'K5272 . ^\ 


"0 


7 


6.82*5 


5 


5tt 
7 


0.7818 


.1 


* 6 


5.*728 


6 


6tt 
7 


q.^339 


3.9628 


•3, 






7 


TT 


0.00Q0 


7.0006 - 


.7 


0 


0.0000 




. 1 ' 








n 


m 00.6990 

• n 




This computation is illustrated, for n 7 , in the last 
column of Table II . 

Remembering that there are n 2 pairs altogether, 
we obtain 

V • 21 

n , n2 . *, 



Exercise 10 

Compute the approximation of ^iven |n (8) : 

t (a) for n = 10 ' 
(b) for n = 100. 

Exercise iN 

(a) Show that the error in estimating p -by (8) is less 

1 n 
* than — . * * 

n 

(b) Show that the estimated value of p^ given by (8) is 
less than the true value*. 

v5 



4.6 Takiqg the Limit 

Let us call the right side of (8) q , We shall show 

tnat £i2 % 9 i"- Since, as you have shown in Exercise 

11(a), |q n - p n | < i, it will follow that lim Pr = |, and 

we shall be done. 

' * « 

"We shall begin with a little algebTa^. Starting from 

the definition of q r 

v " * \ • ~ ■ 

n 



> n sin *— n . 

2Zl = I n 5111 V 



* n n 1 



14 
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Now look again at (3) where we v partitioned the interval 
tO < 9 < tt.. The right hand endpoints of the subintervals 
arfe ir/n, 2*/n, etc. The numbers ~ sin are "the values 
of '.the function f(6) = ^ sin 0 at these endpoinls. And # 

TT fc * 

is the length of each subinterval. In other words 



(see Figure 4) , this is a Riemann sum for the integral 




0 n jr n 2rr . 3tf *i tt 5 tt n 6 tt n 7 tt - 8tt * ' 9 * 

n n r> ~ _ ™ " ' ~ ' 11 — 



Figure k. Relationship between |^-s i n jj- (rectangular areas) 

and r -sin 6 d6 (n » 10). ' v 

Jo tt - 

- sin 0 do. Thus, 

Urn I ICsinlijl-fisi, 

'From (9) we then get 

(10) * lim q n = f 77 i 

• n-*» " Jo 71 



sine d6. 



sin e d6 



i<- c<,s( *'. I. 



Exercise 1£ : 

Suppose In setting up our calculation we measured the angle <j> 
which the needle makes with the northward direction, instead of 8. 
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(a) What inequality would replace (1)? 

(b) Over what values would 4> range? 

a 

(c) What integral would replace the one in (10)? 

(d) Show that this integral equals — . 



5. VARYING THE NEEDLE LENGTH 

We have already remarked that our main, result is 1 
independent of the lengtn L of the needle. But it is 
very jnuch dependent lipon the fact that th^ distance 
between ccrrisecutive lines equals the needle length. 
Certainly if we were to switch to a needle,' say, half 
as long*, while continuing to use the same ruled paper, 
the needle would be less likely^fjo cross a line. 

In^practice it is unrealistic to insist upon this » 

equality. It is very, likely that in setting up this 

experiment you will have available some ruled paper whose 

lines are a distance D apart, and a needle of^length L, 

where L f D. tt turns out that, as long as L is less 

than D, the probability .that the needle will cross a 
2 - T 

1-ine is - x =r. In fact, it f s not at a.ll hard to modify 

. TT u » 

our proof, starting with (1), to get this result. We 1 11 
lea?ve i.t to you (see Exercise" 13) . • * 

Exercise 13 

Suppose the needle has length L but the parallel lines are D 

units apart, where L < D. 

(a) What inequality replaces (0? * 

A (b) What Integral replaces the one In (10).? 

2 L 



Exercise \h 



(c) Show that this integral equals — x — 



The result stated in this section clearly cannot be true if 
2 L 

t > D, since then — x — might be greater than one, and could not 

59 ' . " 



>possIb.!y be the fraction of tosses which are successes. Exactly 
where does our proof break down if we try to modify It as in 
Exercise 12? 



6. GETTING THE COMPUTER TO HELP ' ' 

Whatever the ment of this approach to approximation 
of v, in view of Exercise 8, no one can claim it to be 
a realistic way of getting a good estimate. 'Remember, 
not only does Exercise 8 show that the ^accuracy *after ' 
10,000 tosses cannot be better than three places, byt 
the discussion in- Section £.2 says there is 'no reason " 
to presume it will even be that good. 

It is pos§ible for a computer, figuratively at least, 
to toss the needle for us. In our proof we described the 
approximate position of the needle by an ordered pair 
(y,6) where 0 < y < L and 0 < Q < tt. Since the result 
does not depend on L anyway we .can take L = 1 for con- 
venience. We can ask the computer to pick a number at" 
random between 0 and 1, and call it y; pick another number 
at random between 0 and 4 tan -1 1 (note that 4 tan M l = tt) , 
andtcall it 6; and then determine if 

y + sin e v*i 
is true. ^ * 

Tne machine can easily ^ount the number of "tosses" 
T and of "successes" S$ We can program it to perform 
any predetermined number of "tosses" 'and .then tb compute 
an estimate of tt just as in Section 3.2 / 

In an actual computer run,* we asked the computer/ to 
"toss the needle" 100,000 times. It reported 63,449/ 



* 



Computer work for this unit supported by the University of /Mary- 
land Computer Center. 
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"successes," giving 0.63449 as an estimate of |, an'd 
thus 3.15214/as an' estimate of tt, an error of about % 
of one percent. The BASIC program we used contained 
only 10 linps. If you know BASIC or any other cpmpu^r 
•language, you should be able to try this yourself . ^ ' 

Exercise 15 * * 

I ' * * 

(a) If you know a computer language, write^a computer pcogram 
to simulate Buffon's needle experiment for 100,000 * c I 
"tosses ." 

(b) If you have Recess to a computer, run this progr,am j 
compute the estimate of tt resulting from thWrun," 

^ and compute the percentage error in this' estimate.. 



t 
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7.'. ANSWERS TO EXERCISES 



'2. 
3. 
k. 
5. 
6. 
7. 
,8. 

9. 



10. 

n. 



52 
J_ 
12 

(a) jg (b) 

1 

2 



3* 



(a) 
(a) 
The 



(c).- 



IT — , 3 

17 (b)-^- (c) 0 



best. possible^ result is s - 6366, giving £= 0.636$_c*n.d 
3.1U7. ' .J 



(a) 


j 


1 2 


3 


4 * 5' 


6 


7 


8 9 10* - 




k, 
J" 




1 


0 0 


0 


1 


* 6 10 / 




m io 


=* 68; 


p io s 


0.68. 
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(b) 


J. 


1 2 


3 . 


<*3 




45 


46 ... ,50 




J 


96 t 93 


90 . 


2 


1 


1 


0 0 ' <• 
4 



For 50 < j. < 100, k\ = k ]Q0 

• Soo"" ,2; p,6o *' 

(a) p 1Q 4 0.6314 (b) p 10Q s 0.6366 



j; k ]0Q -= 100; 



You are/looking for an Inequality of the/form 

£nsJn,Ji' 1 I 
' 0_<p n - < n . 

To get It, remember the exact value of d is 
• r n 
n 



1 



n z - I k 
j-1 



Sfftce the formula for Pr involve/ k., write down 
defining Inequality for .k. ; 

k. < n (1 - sin il) < k. + 1. 
* J n j 

Now take, the fol lowing steps. 



the 
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Subtract k. from all, three qualities in the inequality 

0 < n(1 - sin Ji) - k- <\ 
~ n j 

Suln from j — 1 to n, and divide by n 2 - 

n n 
Z n (l - sin Jl) - Z k. 



•0 < M 



^plit the sum of the left, and simplify the far right 

n n 
Z n - Z r 
j«1 j=1 



0 < 



" " . n 

I n - I n sin J £ k : 

n . J 

^i— <I. 

,2 n 



Then split the fraction in the middle this way, 



n - I k. 



J £ n;sin ' , 

0 < LzJ . " < 1 

n 2 n 2 

to get the inequality you seek: 



> Z n s in 



0 < p 



n 



12. 



13. 



14. 



(a) y + L* cos <J> > L. 

(c)K 2 I cos * d+. 

(a) y + L sin e > D. 

* • 

(b) I 1 ]=• sin e d6. 
Jo D 



I- 



The k. defined in (6) is not an accurate count of -those values 
of i from 1 to n for which (5) is falsj. For certain J, ' 
1 " sin n < °'« s0 ^5) is true for all i and false for 0 i . 



But k. < 0. 



63 



15- (a) Here is one BASIC program which should work. Depending 
on tfce computer you are using, it might require minor 
modifications. 

100 RANDOMIZE 

2^ FOR 1 TO 100000 

300 Y'= RMD 

- A00 A « A a RND -fz ATN(l) 

500 IF Y + SIN(A) < 1 THEN 700 

600 S » S + 1 

700 NEXT T" 

~~ "Boo r~=~'s/7o¥o6o 

900 PRINT 'PI = 1 ; 2/R 

1000 END 
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A STRANGE RESULT IN VISUAL PERCEPTION 



1 . INXKODUCTION 

1.1 The Modeling Problem , * 

In 1962 Tom N. Cornsweet reported on the 
experimental verification of some seemingly paradoxical 
results.* He (and others) liad predicted these results 
by constructing a remarkably simple mathematical model 
of the^pyocess by which the eye "sees" bright light, 
based in turn on some ve/y simple physiological assump- 
tions. We -shall describe the experiment and then present 
his model , * 

1.2 The Experiment 

Essentially, subjects whose eyes had had a chanci 
to adapt to darkness* fixated upon a point in a brightly • 
lighted region, across which was a non-opaque bar which 
filtered out a fixed amount of the light. At a certain 
time t 'the bar rather abruptly appeared much brighter 
than the background even though it* was less b 4 rightly„ 
illuminated! Then the apparent brightnesses gradually 
became equal, and final! y^at t*ime t + the background 
became and remained si ightly brighter than the bar. The 
subjects were asked to press buttons at the times t" and 
t as accurately as possible. This experiment was rerun 
for different background brightnesses. T^e details of 
the experiment can.be found in Cornsweet 's article (pages 
261-263), In^figure 1 -average empirical values of t" 
(lower curve) and t* (upper curve) are plotted against a 
measure of background brightness. 



*Too N. Cornsweet, "Changes in the Appearance of Stimuli of Very 
HighauminanceV' Psychological Review 69(1962): 257-273. 



S8 



„* v <-v ju *»u ( OU /O 80 

Brightness (trolands x 10^) 

. Figure I.* 1 Times at which apparent brightnesses of bar and 
background were equal, plotted against background brightness. 
• (Source: Cornswee't, p. 261) * 
i 

2. PHYSIOLOGICAL BACKGROUND 

^ ' - • • «*. 

Before constructing a mathematical model, we must 
have some idea what we are modeling. In this section we 
shall describe briefly what we need to know about the 
physiolo'gical process by which the eye "sees." 

' 2.1 Activation . 

In the rettij^f of one's eye there are certain cells, 
call "receptor" cells. Each of these contains molecules 
which are capable of being "activated 1 ' by light. This 
activation results in the dischaxge of a certain chemical, 
Accumulation of a certain minimal amount of this chemical 
triggers a nerve impulse which results in one "seeing", 
the light . • 

For various reasons^' the chemical tends to be 
destroyed very quickly after its discharge. This means' 
that, if there is to be any hQpe that the required 

Y 
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minimal amount accumulates before it^is destroyed, this 
hope will depend upon very many molecules being activated 
in a shojjf period pf time/ That is, it wall depend on a 
high ram of activation. 



2.2 Regenerat i*an . 

When one 4f these molecules has just been activated 
it is not capable of "immediately" responding to further 
stimulation. It must take a small, but not negligible, 
amount of t,ime to recover, or be "regenerated." So at 
any moment *>niy a certain fraction of the molecules are 
»in this regenerated state.. 

v 3. THE MODEL 



5.1 Assumptions of the Model 

There are* four mathematical assumptions* we must make 
in order to construct our model. We hope each of these 
assumptions wllj seem reasonable to you. As for whether 
they are "true," we cannot say for sure, but we remind 
you that a very striking and ' seemingly paradoxical result 
predicted by the model has been confirmed experimentally. 



(1) 

(2) 
* 

(3) 
(A) 



How much light one "sees" (one's perception) i§ 
directly proportional to the rate of activation. 
The rate of activation is directly proportional 
to the amount of light (brightness), falling upon 
the receptor cells. 

At any time £, the rate of activation is also 
directly proportional to the fraction of the 
molecules which are in the regenerated s^ate at 
that tiine. t 

At any t^me t, activated molecules are' 'being 
regenerated at a ratej directly >proport,ional to 
the fraction of the molecules* whicfi arts v? 
activated Zt that time. * 
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3.2 The Assumptions Rewritten Mathematically 

Let x ■* x(t) stand for the fraction of molecules 
which are regenerated at time t, Then the number of 
regenerated molecules is mx, where m is the total number" 
of molecules. Assumptions (2) and (3) together say that 
an amount of light q shining^ oivthe retina activates 
molecules v (decreasing mx) at the rate cqx, where c is a 
positive constant of proportionality. Assumption (4) says 
that at' the Sams time\other molecules a£e~ being 
regenerated (increasing mx) at the rate k(l-x), where k 
is another positive proportionality constant,. Notice 
that while the rate of activation of regenerated molecules 
depends upon q, the rate of regeneration of active 
molecules does not. , r 

All told, then: * > 

(1) (mx) ■ = kd'-x)" - cqx. 

Finally /assumption (1) says that our perception of 
light is proportional to cqx. Since we are studying ' 
perception, it is this quantity we v are interested in. 



3.3 Summary of the Notat^n v 

A lot of notation is Winning to pile up, and 
there will be more. There is, rio need to remember all 
the notational details. Ju$t keep in mind: 



«q (the amount of light entering the eye) is a posi- 
tive constant controlled by the experimenter* 

c, k, and m are positive constants determined by 
the physiology of the eye I and are not controllable. 

x is a function of the time t, 

cqx is proportional to perception, and is what we are 
Interested in. 
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4. A FORMULA FOR PERCEPTION 

4.1 Solving the Equation of the Model 

<* 

Since -m^ is a constant, Equation (1) can be written 
mx' = k(l-x) - cqx 
mx f r= (-k-cq)x + k 

•(2) X ' = ♦ £. 

m m 

We are getting buried by notation. At the price of 
introducing still more letters, let's simplify it. In - 
-effect, Equation (2) just says x' is a linear function 

(3) x' = rx + s \ 1 

where r and s are constants. . 

4 Here is a very useful tricjc for solving any equa^g^n 
Jike Equation (3). You may want to remember it. Just 
make the substitution * 

z = rx + s . %m . 

Then z' * r*' = r(*x+s) 4 = rz. And of course the solution' 
of z 1 = r2 is 2 * Ce rt . .Going back to the x -notation, 
this is rx + s * Ce rt ; or 

(4) ' = C*e rt - i 
whe^e C* * . " . 

" 4,2 Determining the Constant * 

To determine C* we must. know the value of x for any 
one particular t. Remember, in the experiment the sub- , * 
ject's eye was allowed* to adapt to darkness before being 
exposed to a bright light,* "In other words^, almost all 
the molecules in the receptor cells were in the regenerated 
state just before the initial exposure, so at that moment 
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x = 1. If we label the moment of first exposure t = 0, 



then x(0) * 1. 
Q 



Putting t = 0 an<#x * 1 in Equation (4) yields 



C* * 1 * i * r + s 



r ■ r 
and so 

r > s 'Jtt- 



-trr 



4.3 Determining the Activation Rate 

To interpret Equation (5) we must return to our 
origina} notation. Comparing Equation (3) with 
Equation (2) we see 'that - 

* 

r = - k - c q 



k 

and s - so 
m 



x cq -(k+cq)t/m k 

X F+cq e . < ? k + cq 



and the rate of activation is 
(6) , cqx 



2 2 f 



„ k + cq k + C q * 

The right side, of Equation (6) tells us, if our model 
has any validity, how bright a light will appear to one 
as a function of time t, assuming it is turned on at time' 
t = 0, and that one's eye^was in darkness before then. 
Let's give this function a name: f(t), and let's 
investigate -its properties. 
* • * i * 

. ' 5. ANALYSIS OF fHE FORMULA 

5.1 Simplifying th,e Notation 
r 

Again, let's simplify the notation to avoid being* 
buried by it. Let's write 
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M = X^—> N = > k * cc i 
k + cq* m 



Then we have 



P = 

k + cq 



(7) 



f(t) - Me* Nt + P. 



There vis no point in trying to interpret* M, N, and 
P physiologically. What we should remember is the 
material in the box in Section 3,3, from which it follows' 
that M, N, and P are positive constants,^ 1 

5.2 The Shape of the Perception <€raph 

It is easy to confirm (Exercise 1) that f(0) = M + P, 
lim f(t) = P, f'(t) < 0 for all t, and f"(t) * 0 for all t, 

t-»® -j 

Therefore the graph of f(t) for positive t looks like 
Figure 2. The horizontal asymptote at P represents an 
equilibrium position. It corresponds to that value of x 
(the fraction of molecules which are regenerated) for 
which regenerated molecules are being activated and 
activated molecules are being regenerated at the same rate. 




Figure 2. A typical graph of the activation rate as a 
function of time.. 
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5.3 The Effec/'Of Changing the Brightness 



ec/'Of Chan* 
ftange q, th< 



If*we d&ange q, then we a^lso change M, N and P. But^ 
they are still positive ^constants* So w,e get a new f(t^*. 
whose graph has the same shape as that, in Figure 2, but 
with* a different intercept and a different horizontal, .* » 
asymptote. ( ^ 

7 

6. RELATION OF THE FORMULA TO THE EXPERIMENT , 
6.1 Explaining the Experimental Result 

v Now let's look at the experiment. What happened 4 
was that the subject's eVe received two amounts' of light i 
q^ from the background, ajid q2 < q^ through the bar. So 
there were two graphs like the one in Figure 2: 9 



f A (t) = Mje^l^ + P 



tfith intercept M^ + P^ and asymptote at Pj, representing 
the subject's perception of the background, and 

f 2 (t) = M 2 e" N 2t - m ^ 

with intercept M 2 + P 2 and asymptote at'P 2> representing 
the subject's perception of the bar. 

In Figure 3 we have drawn two curves of the right 
shape on tjie same pair of coordinate axes. We have not 
considered the specif rc values of the' ^constants . We % 
just want to show that it is at least believable that 
the two curves might intersect twice in the region t > 0», 4 
in which- case there would be two reversals as to which 
cui»ve was the higher of the two. If this actually hap-" 
pens', then for t" < -t < t the subject will perceive the. 
less brightly l.lluminated bar to be brighter than the 
background, and there will be^ two reversals of apparent Va 
relative brightness of bar and background. Further, 
'Figure 3 suggests thar the first reversal (at t = t ) will 
be rather abrupt ," whi le the second (at t - t + ) will *be 
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much more gradual. This is just what happened in the 
experiment ! 

. f(0 




Figure 3. Two typical graphs of the activation rate (same 
subject, different brightnesses). 

iL? — Predicting the Experimental Result 

Now for that' all-important -"If Could we have 

actually calculated that f^t) and f 2 ( t ) intersect ' (and 

thus have predicted the experimental result), rather 

than settle for "it is at least believable . /> .and 

(after jhe fact) "This is . . . what happened'^ 

« 

, ■ 'This- amounts to solving f^t) = f 2 (t), or 



(8) 



l^e-Nl* 



=.M 2 e" N 2 t t p 



for t. -There is usually no easy way to do this. But 
for certain carefuWy selected values of q 1 and q 2 it 
can be done. (Remember, q is the only thing We can 
control). Let's select qj and q 2 so that -Nj = 2N 2> . 
(This can be done. In Exercise 7 you are asked to show 
^that qi = Rq^ will do it, where 
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k + 2cq 2 , 

* = — + 2.) 

cq 2 cq 2 > 



Then Equation (8) becomes 



M 1 e" 2N 2 t + P 1 =' Nke'V + P„ 



(9) M 1 e" 2N 2 t - M 2 e" N 2 t + (P - P ) = 0. 



Putting y = e" N 2 t we get 



(10) M x y 2 - M 2 y + (? l - P 2 ) - 0 

which is a plain old quadratic equation. (This is the 
payoff for making ■ 2N 2 .) Its solutions are of course 



(11) y = e V ^= -i 2 w , M 1 1 IL. 

1 . _ 

There will be two positive value! of- t satisfying 
Equation (9) if there are two values of ybetween 0 and 
1 satisfying Equation (10). So we must show that the 
right side of Equation (11) is (a) real, (b) greater than 
zero, and (c) less than one. 

The sticky part of this is (a). Some fairly messy 
algefira is required, which we have put in the Appendix. 
It turns out that the right-side of Equation °(U) is ' • 
real if q 2 is lar§e enough compared with k/c . After 
that, (b) and (c) are relatively easy, and we ! ll leave 
them as exercises (see Exercises 5 and 6). 

So if q 2 is sufficiently large (see the Appendix), 
and if q 1 is suf f iciently 'larger than q 2 (Exercise 7), 
then ^(t) and f 2 (t),will intersect at two positive 
values of t, producing the paradoxical effect confirmed 
by the experiment. 
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APPENDIX 



This appendix is devoted to showing that the 
solutions to Equation (10) are real' provided q 2 is large 
enough. As stated in Section 6.2, these solutions are 

M 2 t A\ 2 l - 4M 1 (P 1 -P 2 ) 

y = ' 
•« * 

They will be real (and distinct) i£ 

» - 

M 2 2 - 4M 1 (P 1 -P 2 ) > 0. 

In working with this inequality, we shall use 

\ 

(a) the definitions: j 



' 2 2 
C q 1 

1 k+cq n 



(b 



2 2 
c q 2 

K + cq 2 



cqjk 



P l" = F^cqJ 



cq 2 k 
k+cq^ ' 



the fact that q^ = Rq 2 where 



R = 



k+2cq- 



cq 2 cq 2 



♦ 2; 



(c) the. fact that k + cRq 2 = 2(k+cq 2 ), which you are 
asked, to prove as an exercise (see Exercise 4). 



We have 



tby (a)) 
(by (b)) 
(by (c)) 



- 4Mj (P 1 
















' cq^k CQ Li^ ' 






(k+cq 2 ) 2 


k + cqj 


k*cqj k+cq 2 






4 4 . 
c q 2 


4c 2 R 2 q 2 2 


' cRq 2 k cq 2 k 




(k + cq 2 ) 2 


k+cRq 2 


k+cRq 2 k+cq 2 




c q 2 • 


2c 2 R 2 q 2 


2 


cRq 2 k 


cc» 2 k' 


(k + cq 2 ) 2 


k + cq 2 


) 


2(k + cc( 2 ) k + cq 2 



0 
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4 4 



c q 



,c 3 R 3 q 2 3 k + 2c 3 R 2 ^q 2 3 k 



. f ' (k+cq 2 r 

c 3 q 3 

= ^(cq 2 - + 2R 2 k] . 

V . 

For this expression to be positive, the part in 
parentheses^must be positive, ing ttfe definition of R 
in (b): 

4 

cq 2 - R 3 k + 2R 2 k 



k+2cq 2 


3. 


k+2cq 2 


2 


k + 2 

( 


k 


[ cc *2 J 


I Cq 2 



cq 2 



; TT 4 + 6k 3 cq 2 + 12k 2 c 2 q 2 2 + 8kc 3 q 2 3 2k 3 + 8k 2 cq 2 + 8kc 2 q 2 2 

c 4 q 2 4 -k 4 -6k 3 cq 2 -12k 2 c 2 q 2 2 -8kc 3 q 2 3 +2k 3 cq 2 +8k 2 c 2 q 2 2 +8kc 3 q" 3 

C q 2 . / 



c 4 q 2 4 .k 4 -4k 3 cq 2 -4k 2 c 2 q 2 2 % 

3 Z r~ * 
, C q 2 

For this fraction to be positive, the numerator must be 
positive. That is, ^ * 

A « c 4 q 2 4 - k 4 - 4k 3 cq 2 - 4k 2 c 2 q 2 2 > 0. 

Since it is q 2 which we can control, let's rearrange 
terms : 

t 

c 4 q 2 4 . - 4k 2 c 2 q 2 2 - 4k 3 cq 2 > k 4 . 

The left side of this inequality's a fourth degree v poly- 
n,omial in q^. Since th«e leading coefficient c 4 .is 
positive, the left side approaches +« as q 2 approaches 
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+*. And so the. left side 'must be greater than k 4 for 
q 2 large enough. 

Actually, you can easily check that if q 2 ' = 3k/c, 
then the left side equals 33k 4 , which is certainly greater 
than k 4 . 



EXERCISES 

1. Starting with. Equation (7) : r ^ 

(a) compute f'( t ) and f"(t); 

(b) show that f'(t) < 0 for all t; 

(c) show that f"(t) > 0 for all t; « 

(d) .show that f(0) » M + P; 

(e) show that lim f(t) - P. - 

2. Find the time t between t" and t + on Figure 3 at which the 
second light is perceived to exceed the first light 'in brightness 
by the greatest amount. 

3. In Equation (2), find the value of x for which we have equilibrium, 
and show that this is consistent with the remarks in Section 5.2' 
about the equilibrium position. 

f 

4. Show that if " 

k + 2cq 

R - £ = JL. + 2> 

cq 2 cq 2 

then 

k + cRq 2 - 2(k + cq 2 ). 

5. Show that the right side of Equation (11) is positive,- assuming 
it is real* 

6. Show that the right side of Equation (11) is less than one, 
assuming it is real. 

'• /« 

/ 
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ANSWERS TO EXE RCISES 



(a) 

* 


f'(t) » 


(b) 


e z > 0 


(c) 


same as 


.(d) 


f (0) - ] 






(e) 


lim tMe" 








since e 



•Nt 



-Nt x 



P = 0 



Nt 



-* 0. 



log - log M 2 N 2 

/>-"> . - •• 

(Maximize* the funptW g(t^. = f 2 (t) - f^'t) by setting 
g'(t) =» 0 and solving for t. ) 

X S k + cq ' (Set X, * = 0 and . * So * ve for x -> 

For this x, the rate of activation woujjd be 
cqk 



cqx 



= P,. 



k * cq 

as stated in Section 5*2. 



k + cq. 



k + cRq 



k + cq 2 . k + cq 2 * 



}c + »2cq 2 



k + *|k + 2^- v 
k + cq. 



- 2. 



The solution with the plus sign is certainly positive, 
other will be positive'if 



\ " ^2 2 " 4M l (P l - V.* °' 



^ > ^ 2 2 - AM 1 "(P 1 - P 2 ), 
m 2 2 > - 4fl 1 (P 1 r p 2 ). ; 
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For tMs^o be true we require P 1 - P 2 > 0. To see that this 
is so , we write 



cq,k cq,k 
p p _ ' 1 



1 2 k+cq 1 k+cq 2 } 



2 2*2 2 
k cq x +c q x q 2 k-k rq"-c q^k g 



(k+cq 1 )(k+cq 2 ) v 

k^c ( . 

(k+cqJCk+cqO ' IW » 

.y. ' ■ - 

which is positive since q > q 0 . 
V 

We^aist show 



, M, 2 - ^(P^) < 4^ - ^ + M 2 2 

- ' P 'l. : P 2 .> "2 " "l' . 

+ M 2 ^^P 2 . ' 

To show this, note that ' ' • - 

,Jf , 2 » 2^ 



f". " w„, -fair 



« *c q, +cq,k 

ii 

and similarly that 



Remember that- q > q^.v 
1 . t 



•1 



J.. 



7.» N x k+cq x <.k+cRq 2 » * 

" ~ " 2 (us-rft^K^ercise 4)\ 



/ 
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N 2 - k^ k+cq 2> ) U ™'* ^" „ . - *" 
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